series mu=b*p/(q-1)
series gini=2*@beta((2*p),(2*q-1))/(p*(@beta(p,q))^2)
series theil=log(mu/b)+@digamma(q)-@digamma(p)

